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The model-independent "box" parameterization of neutrino oscillations is examined. The invari- 
ant boxes are the classical amplitudes of the individual oscillating terms. Being observables, the 
boxes are independent of the choice of parameterization of the mixing matrix. Emphasis is placed 
on the relations among the box parameters due to mixing-matrix unitarity, and on the reduction of 
the number of boxes to the minimum basis set. Using the box algebra, we show that CP-violation 
may be inferred from measurements of neutrino flavor mixing even when the oscillatory factors have 
averaged. General analyses of neutrino oscillations among n > 3 flavors can readily determine the 
boxes, which can then be manipulated to yield magnitudes of mixing matrix elements. 
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• I. INTRODUCTION 

, If neutrinos have mass and are non-degenerate, then their flavors may oscillate as they propagate. Resonant 
^ ^ oscillations for the sun |l[ , oscillations for the atmosphere [|| , and the LSND data each require a different neutrino 
mass-squared difference if neutrino oscillations are to account for all features of the data Since three-neutrino 
1 models can have at most two independent mass-squared differences, a sterile neutrino is apparently needed to reconcile 
all the data while retaining consistency with LEP measurements of Z — > i/F Q . Several four- neutrino analyses appear 
I , in the literature |^J^. It is also possible that some data will turn out to have an explanation other than neutrino 
^ ■ oscillations, in which case three-neutrino oscillations may be sufficient. So our task is to examine the physics of 

, neutrino oscillations with three or more mixed flavors. 
0^ ■ Oscillation probabilities depend on products of four mixing-matrix elements. Several parameterizations of the 
\l mixing matrix in terms of rotation angles have been introduced, beginning with the pioneering work of Kobayashi 
and Maskawa [Q. With three or more neutrino generations, the oscillation probabilities are complicated functions of 
the neutrino mixing angles. But oscillations are observable and therefore parameterization-invariant. One must ask 
qI^ , if there is not a better description of oscillations which avoids the arbitrariness of angular-parameterization schemes. 
Recently, we introduced a "box" parameterization of neutrino mixing valid for any number of neutrino generations 
1^ . Oscillation probabilities are linear in the boxes, enabling a straighforward description of oscillation data. Here we 
II ■ present the algebra of the boxes and the unitarity constraints on that algebra. Then we illustrate the boxes' reduction 
Qh! to a basis in the case of three generations, thereby setting the framework for a future phenomenological analysis. 
^ ■ The probability for a neutrino to oscillate from Va to Vfj is given by the square of the transition amplitude: 
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where n is the number of neutrino generations. 



= ^ (0» - 0j) = ^ {E'lU - PiX^ - Ejtj + pjXj) , (2) 

and Vai is the mixing-matrix element which connects the a^^ charged lepton mass eigenstate and the i**^ neutrino 
mass eigenstate. For relativistic neutrinos, is given by 

^ij « ^ X, where Am^ =mf — m?j. (3) 
With a little bit of algebra, the oscillation probability may be brought into the form 

P (x) =-2j2J2 MVc.^Vp*^:^ Vpj) sin^ (4) 



i j^i 

J2 E MVa^Vp\V:^ Vf3,) Sin (2$,,) -f 5^p. 
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The probability for an antineutrino to oscillate from Va to Vfj is obtained by replacing V with V* . This is equivalent 
to changing the sign of or the second term in equation (^. 

With the familiar case of two neutrino flavors, the mixing matrix V has the simple form of a rotation matrix (phases 
cancel in oscillation probabilities for Majorana neutrinos, and may be absorbed into the definitions of Dirac fcrmion 
fields): 

Tr_f COS 6 —am6\ , , 

^ = I sine cose ) ■ 



The oscillation probability in the two-flavor case is simply 



^12 



P (x) = SaB + sin^ 29 sin^ —r^x , n = 2. (6) 
- \ J 

The mixing-angle parameterization is a natural choice in the two-flavor situation. 

The formalism becomes more complicated with three flavors. An arbitrary 3x3 unitary matrix has three real 
degrees of freedom and six phases, but 2n — 1 = 5 phases may be absorbed into field redefinitions. The original choice 
of the four remaining parameters, due Kobayashi and Maskawa to describe quark mixing, is |Q 

S1C3 S1S3 
C1C2C3 - 52536*'' C1C253 -I- 52036"^ I , (7) 
C152C3 + 02536**' C1S2S3 - 02636"^ 

where Ca = cos 9a, and Sa = sin^^. There is arbitrariness associated with the placement of the phase, since we absorb 
five relative phases into the field definitions. Because of this arbitrariness, the phases of individual matrix elements 
are not observable. 

The observable oscillation probabilities are quite complicated functions of the angle-based parameterizations. As 
an example, consider the product V22V'2*3y3*2V33 appearing in the j/^ — )■ i^r oscillation probability: 




V22V2%V*^V33 = clsl [slcjist + Gcj + 2c\ cos 25) - c 
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^{l + cl){cl~sl){sl-cl)coi5 + iJ, n = 3. (8) 



s 



where the Jarlskog invariant J |9j has the form = 015^02520353 sin 5 in this parameterization. 

The expression (but not its value) on the right-hand side of equation (^) is convention-dependent, as well as being 
unwieldy. Our development of a model-independent parameterization is motivated by the arbitrariness and complexity 
of this traditional approach. 



II. THE BOX PARAMETERIZATION 



The immeasurability of the individual complex mixing-matrix elements in the quark sector has been addressed by 
numerous authors |p|-p^. Measurable quantities include only the magnitudes of mixing matrix elements, the products 
of four mixing-matrix elements appearing in the oscillation probabilities, and particular higher-order functions of 
mixing-matrix elements 11,1^]. As evidenced in equations (|l|) and (^), neutrino oscillation probabilities depend 
linearly on the fourth-order objects, 

^ V^,V^pVp,Vl = V^,V:^V;,Vp„ (9) 

which we call "boxes" since each contains as factors the corners of a submatrix, or "box," of the mixing matrix. For 
example, the upper left 2x2 submatrix elements produce the box 
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□ 22 = Fiit/i*2 1^2*1^22. (10) 



The name "box" also seems appropriate in light of the Feynman box-diagram which describes the oscillation process. 
Examination of equation (^) reveals a few symmetries in the indexing: 

-□/3, = = P'ul^ - "^n^, . (11) 
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If the order of either set of indices is reversed {id est, j ^ i or f3 ^ a), the box turns into its complex conjugate; 
if both sets of indices are reversed, the box returns to its original value |0. And if V is replaced by , then 



trn* 



Boxes with a = f3 or i = j, are real, given from equation (H) as 
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— \Vai\'^\Vf3i 



and 



(12) 



We call boxes with one and two repeated indices "singly-degenerate" and "doubly-degenerate," respectively. Boxes 
with a ^ P and i ^ j are called "nondegenerate" . As can be seen from equation (|4|), singly-degenerate boxes with 
repeated flavor indices enter into the formulae for flavor-conserving survival probabilities, but not for flavor-changing 
transition probabilities. Degenerate boxes with repeated mass indices (including the doubly-degenerate boxes) do not 
appear in any oscillation formula. Degenerate boxes may be expressed in terms of the nondegenerate boxes, as will 
be shown shortly. This possibility and the symmetries expressed in equation ( pT| ) allow us to express combinations of 
boxes in terms of only the nondegenerate "ordered" boxes for which a < (3 and i < j. 

Using the symmetries expressed in equation (11), the oscillation probabilities (^ in terms of boxes become 
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'ifSj sin2<I>jj] 



(13) 



where we have defined the shorthand 



(33 



Re I 



P3 



and 



Im I 



Pi 



From equation (11) we 



deduce that the Js are antisymmetric in both flavor indices and mass indices; Rs are symmetric in both. Survival 
probabilities P {x) = 1 — "^g^^P (x) are more simply expressed in terms of degenerate boxes, or |l^|s, rather 



than nondegenerate boxes. From equations ( p^ ) and (12), they are 
P (x) = 1 



4EE 

i—l j>i 



sin2 = 1-4^^|F„, 

1=1 j>i 



|14,f sin2 



(14) 



Interchanging a ^ P in equation ([13) gives the time-reversed reactions P (x) : 

1^0 I'a 



P {x) = (5„/3 - 2 ^ ^ [2 "'R^, sin2 + J^, sin 2$,,] 

■i—l jyi 



(15) 



Ignoring possible CP-violating phases in the mixing matrix, the number of real parameters determining V is the 
number of rotational planes available in n~dimensions, N = \n{n - 1) . Deter mining these N parameters determines 
the complete mixing matrix. Conveniently, there are N transition probabilities P (x) — P (x) . Thus, all of the 

I^Q -t 1^13 Vf, -* Va 

information in the mixing matrix is contained in the N transition probabilities. In this sense, they form a convenient 
basis for determining all oscillation parameters. Of course, if the same transition probability is measured at two or 
more different distances, then all N transition probabilities may not be needed to determine V . 

Allowing CP-violation in the mixing matrix, there are N real parameters and \{n — l){n — 2) phases, for a total 
of (n — 1)^ parameters. With CP-violation, however, there are 27V = n{n — 1) independent transition probabilities 
P (x) . The number of transition probabilities exceeds the number of independent parameters, so they again form a 

convenient basis for determining the mixing matrix. In reality, only the three flavor indices e, /i, r are easily accessible. 
Moreover, some of the N parameters in the mixing matrix, namely those which rotate sterile states for n > 5, are not 
accessible at all, which complicates the counting. 

The transition probabilities for which a 7^ /3 in equation ( |l3|) may be conveniently expressed in matrix form. The 
matrix of boxes is an A'^ x TV matrix. For three flavors, we have 



Tin = 3) = 



where 
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-4 Re(S) S'2($) + 2 Im(S) S'(2$), 



sin $12 

and 5''($) = ( sin*"' $23 
sin*" $13 



(16) 



(17) 
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For the time-reversed channels, or for the antineutrino channels, the sign of the lm{B) term is reversed. The box 
parameterization is especially well-suited for considering higher numbers of generations. The matrix B merely acquires 
extra columns when new flavors are introduced; extra rows are not accessible at energies below new charged-lepton 
thresholds. Furthermore, oscillation probabilities are linear in boxes, no matter how many generations. 

Neutrino oscillation experiments will directly measure the boxes in equation ( p^ ) , not the individual mixing matrix 
elements, Vai- But one would like to obtain the fundamental Vai from the measured boxes. We develop here an 
algebra relating boxes and mixing matrix elements. 

Some tautologous relationships between the degenerate and nondegenerate boxes are easily confirmed using equation 
(|); they hold for any number of generations: 



|2n/ |2 _ ain _ 



= »n^, = - , and (19) 



(rj^a^ 13, andx^iy^j). (20) 



Equations ( |l8| ) and (|l9| ) are themselves special cases of the more general 

-'^os, = [v^.v:^vp,v^,] [v^,v;^vs,v;^ (21) 

and the analogous relation "'D/jj "'^□/s; — "^D/s/ "'^^(ij ■ The relations above hold for both degenerate boxes and 
nondegenerate boxes. 

Due to the symmetry "'D/Sj — > '"l^j/j when ^ tj^ere will generally be analogous but distinct pairing of our 
box equations, differing only in whether the degeneracy or sum is over a flavor index or a mass index. In the following 
we will mainly show only one equation per analogous pair, for reasons of space limitations in this proceeding. 

We may express \Vai\ = ( "'Dai ) " in terms of three singly-degenerate boxes by setting a — f3 in equation (|l9|). 
Then, using equation ( p^ ) to substitute for the singly-degenerate boxes yields an expression for the doubly-degenerate 
box in terms of nine nondegenerate boxes: 

ai] — I ai] — i ct^rn ^^[""1 ^vn r~i 

ly |4 _ ai|-| _ ^ax ^ay _ '-'ti '-'ax '-'pi '-'Cv ^ f^V ^22) 

^ay ^ax '-'(y "'-'ll}x ^P-y 

where the index constraints are r^a^a, C.=/=p=/=a, fiy^LU=/=a, and x =/= y ^ i. In the three-generation case, 
equation ( ^2|) is uniquely specified by the index constraints. For example, 

IT/ |4 _ '-'■22 I-I23 '-'33 '-'32 '-'33 

- -□*3 12033 -^33 -^^2 ^ ^ 

holds with any number of generations, but it is the unique 5 on 4 box representation of |Vii|^ in three generations. 

We note that all of the relationships in this section follow from the definitions of the boxes in equation (|^) and so 
are valid for any matrix, unitary or otherwise. The constraints of unitarity will provide us with expressions for 



which are easier to manage than the expression in (22) above. 



III. UNITARITY RELATIONS AMONG THE BOXES 

Unitarity requires that 

n n 

T.^v^VnJ^^^l^ and J2v^yVp*y^S^f,. (24) 

j)=l y=l 

Multiplying the first equation in ( p4[ ) by V^^Vxj and the second by V*^Vf3x gives the unitarity constraints for the 
boxes: 
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J2 "'Daj - V^'OmS^,, (25) 

ri=l 



and its analogue. Isolating the manifestly degenerate boxes from the nondegenerate boxes, equation (Pq) becomes 



J2 "'Da, = ^^^^6^J - ^'Daj . (26) 



Summing equation (p5|) over A in the i ^ j case, we find 



o = EE "'^^^ +2EE "'R^^- (27) 

A=lj)=l A=l A=l77<A 

The double sum is over Rs only, since the first term is manifestly real. The resulting conditions on the Js are found 
in equation (|2|) below. Comparison of equation ( |27| ) with equations ( |l7| ) and (|l^) reveals an interesting property of 
the matrix B: 

n 

^ Re{B) = ~-J2\V,,\'\V,,\', (28) 
column o/B A=l 

where the sum is over a column of B specified by fixed i and j. There is an analogue relation for the sum over a row 
of S. 

The unitarity constraint ( p5| ) holds independently for the real and imaginary parts of the sum. We will first 
explore the implications of these constraints for the imaginary parts of boxes, before turning to the more complicated 
constraints for the real parts. The right-hand side of ( P5[) is manifestly real, so the imaginary constraints are simply 

J2 "'Ja, = 0, and "'J/5. = 0. (29) 

Equation (LI) indicates that is an antisymmetric matrix in the indices rj and A for fixed i and j, and vice versa. 

Equation (2£) shows that the sum of elements along any row or column of that antisymmetric matrix equals zero, 
whether the sum is over mass indices or flavor indices. 

Summing the first equation in ( p9|) over A gives zero trivially since a sum of all elements of an antisymmetric matrix 
vanishes by definition. Hence, for fixed (i, j), the first equation in ( p9| ) expresses n — 1 constraints. Thus, the number 
of independent flavor pairs on Js after implementing the constraints of equation ( p^ ) is — (n — 1) = i (n — 1) (n — 2) . 
Ditto for independent mass pairs, so the number of independent Js after implementing both sets of constraints is the 
product j{n — 1)^ (n — 2)^. 

In three generations, this number of independent Js is one. Each sum in equation ( p9| ) contains only two terms, 
leading to 




Im(B) = J J -J \, n = 3, (30) 



with J = ^^J22 0. One consequence of the equality of all |J|s in three generations is that if any one is zero, 
then all "*J/3j vanish and there can be no CP-violation. 

We now consider the real parts of the constraint (|2^), focusing first on the homogeneous constraint for which the 
Kronecker delta is zero. This constraint gives the singly-degenerate boxes as sums of ordered boxes: 

\yx^\^\y>.,\^ = "'aA, = - E - E "'i^^^' * ^ ^- (^i) 

This linear relation complements the relation expressed in equation (|l^). For three generations, each of the sums 
contains two terms, allowing us to express the singly-degenerate boxes in terms of two nondegenerate boxes which are 
measurable in neutrino appearance oscillation experiments. 

The real unitarity constraint (|3l|) greatly simplifies our expressions for a doubly-degenerate box "'Dai — \ Vai\'^: 
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(32) 



where the first equahty is due to equation ([iq ) with j — i. Applying equation (^) to three generations, one finds that 
doubly-degenerate boxes are expressible in terms of the real parts of six ordered boxes, rather than the nine complex 
boxes used in equation For example. 



22 
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^R33 



n = 3. 



In cyclic coordinates a,/3,7, and with x ^ y ^i, 



■yyj 



(33) 



(34) 



When considering n > 3, each sum has more terms, but all terms in the numerator in equation ( |32[) always contain 
Rs to the second order, while the denominator terms contain only the first order of Rs. Thus these expressions will 
be much more manageable than equation ( p2|) which exhibits the fifth order of complex boxes in the numerator and 
the fourth order in the denominator. 

For fixed (j,j) in equation (^l|), A can take n possible values, implying n constraint equations. N ordered nonde- 
generate boxes appear in these n equations. Thus, for N < n, which is true for n < 3, the unitarity constraint ( |3l| ) 
may be inverted to find a nondegenerate box in terms of singly-degenerate boxes. Manipulation of equation (|3^) gives 
an expression in term of the flavor triad (a, /3, 7): 



(35) 



It is known that knowledge of four \V\s, completely specifies the three-generation mixing matrix, provided no more 
than two W_\s are taken from the same row or same column ]l5|. Here, we can use three-generation unitarity and 
equation ( p5| ) to re-write "'R;3j in terms of just four \V\s. The result is 



R/3j - o 



(36) 



which for n — 2> expresses the real part of the box Re 



in terms of the magnitudes of the four complex 



Vs which define the box. Three-generation unitarity may be used again to replace the first five terms on the right-hand 
side of equation ( ^ ) with — |V^fe|^. 

Summing equation ( ^l|) over j ^ i yields another expression for in terms of nondegenerate boxes, which 

further complements equations (p2|) and (p2|): 



The explicit solution of this equation, valid for any number of generations, is, 



(37) 



1± J1 + 4EE "'R^^- 
V i^^ n^>^ 



(38) 



which yields \V\i\'^ in terms of {n — 1)^ Rs, but subject to a two-fold ambiguity. 

We may use the real homogeneous unitarity condition (|3^) along with the tautology (|l^) to obtain constraints 
between nondegenerate boxes, thereby reducing the number of real degrees of freedom. Substituting the tautology 
(O) into the unitarity constraint (^ll) gives 



'—'aj '— ' 



(39) 



This unitarity constraint interrelates imaginary and real parts of n different boxes for any number of generations. For 
example, taking the imaginary part of equation (^) leads to 
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Taking the real part of equation (^) leads to 

"'R^, "'Ra, + "'Raj "I^":' = "'J'-J "'J^J' ^ ^ A ^ a, i ^ (41) 

One may also use the pairs of equations ( |40| ) and (^]) to eliminate the sums and isolate a single R: 



i T air? Pij) i ai t t 8i-\ 

air, _ J;3j RAj RAj + J/3j JAj ^ J \j , .„^ 

"Oa, '''Ra, - -Ra, ^^Jx, ' ^ ^ 

with f3^X^a, i^j. Input from these unitarity relations among Rs and Js is necessary to estabhsh the minimum 
set of independent box parameters. 



IV. INDIRECT MEASUREMENT OF CP- VIOLATION 

Suppose CP is conserved. Then "'^Jpj = for aU index choices. The inference from equation ( ^l|) is that 

"^R^, "'Raj + "^Raj Y "I^-J =0, (^; ^ A ^ a, and z ^ j). (43) 

If this relation is violated, then so is CP. 

For three generations, ^^Jaj = "'Jaj by equation (^), and equation ( ^ may be solved for J"^ directly: 

J2 = "»R^^. /5»R^^. + -R^^. "^R;,^. + -»R^^. /3»R;,^., „ = 3. (44) 

Equation ( ^ ) says that the three real elements in any row (or any column in the analogue equation) of the matrix 
B may be summed in their three pairwise products to yield the CP- violating invariant J^. These real elements on 
the right-hand side of this equation are measurable with CP-conserving averaged neutrino oscillations. Thus, even 
if CP violating asymmetries are not directly observable in an experiment, the effects of CP violation may be seen 
through the relationships among the real parts of different boxes, which are determinable from averaged flavor-mixing 
measurements! Note that if CP is conserved and J is zero, then equation ( ^ ) also tells us that all three Rs in any 
row (or column) cannot have the same sign. 



V. INHOMOGENEOUS UNITARITY CONSTRAINTS AND A BOX-BASIS 

The inhomogeneous unitarity constraints with the Kronecker delta nonzero in equation ( p5| ) are necessary to provide 
the desired normalization of the Vai or the boxes. The inhomogeneous constraints are functions of degenerate boxes 
and therefore purely real: 

+ J2 "'^v^ = (45) 

This equation can be rewritten strictly in terms of nondegenerate boxes by using the homogeneous unitarity constraints 
( ^ ) to replace the singly-degenerate boxes, and equation (|2|) to replace the doubly-degenerate box: 

SaS, + v/-SAS^Sr + SrSr,. = , (46) 

with Sa = Ea^o "Ra., S„ = Ea^a "'R-y; = Er^a ""R-y; = E^^a E.^i "'R^., and 2; ^ y ^ i.^ These 
inhomogeneous unitarity constraints do not involve the Js. Isolating the square root and squaring the equation, we 
get polynomial equations of degrees three and four in the Rs, each relating n(ri — 1) Rs. 

We provide here an example of a basis construction, obtained by substituting in the unitarity equations derived 
above. The unitarity constraints among the Js, given in equation (p9[), are linear and therefore the simplest to 
implement. These constraints may be used first to reduce the number of independent Js to j{n — 1)^ (n — 2)^ . 
Further reduction to independent Js and Rs requires the nonlinear constraints. The homogeneous constraints ( |40| ) 
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and ( p3| ) are much simpler than the inhomogeneous constraints (|6|), but the inhomogeneous constraints must be 
invoked at least once (Otherwise, the the boxes and the matrix element magnitudes \V\ could not be normalized.). 

For three generations, one begins with nine Rs and one J, and seeks a basis of just four elements. Rearranging 
the three-generation equation (|4J) yields expressions for one R in terms of two other Rs and J7. This equation may 
be used three times to eliminate ^^R23, ^^R32 and ^^Rss. The utility of the analogue equation is exhausted to 
eliminate ^^Rss and ^^Ras. As expected, one must next turn to the inhomogeneous constraints ( ^6|) to eliminate 
the last degree of freedom. We are left with a constraint which is quartic in all five of its parameters A = ^^R22, 
B = 11R23, C = 21R32, D = 22R33, and J^: 

= + Bf {A + Cf {BC + BD-CD + J^f 

+ {AD + BD-AB + J'^f [{A + B){C + D) + C^ + j'^f (47) 
+ {A + B){A + C) {BC + BD-CD + J^) {AD + BD-AB + J^) 
x[C + D + 2{{A + B){C + D) + C'' + J^)] . 

We may eliminate any one parameter by either algebraic or numerical means, leaving us with the desired four 
parameters as the basis. 



VI. SUMMARY 



Neutrino physics has entered a golden age of research. New experiments all over the globe promise an unequaled 
amount of data from the sun, the atmosphere, accelerators, supernovae, and other cosmic sources. The latest data 
suggests that more than three neutrino flavors may participate in neutrino oscillations Q. Analyzing such refined 
data requires a consistent, model-independent approach which may be easily applied, and easily extended to higher 
generations. Here we have discussed such an approach, wherein one works directly with the observable coefficients 
of the oscillating terms. From unitarity of the mixing matrix, we derived relations among these CP-conserving and 
CP-violating coefficients for the various oscillation channels. One result which we view as particularly noteworthy is 
that high-statistics data on averaged oscillations are sufficient to determine the conservation or non-conservation of 
CP in the lepton mixing matrix. This indirect test of CP can be traced back to unitarity of the mixing matrix, but 
in the present formulation there is no need to even mention the mixing matrix. 
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Physics, under Grant No. DE-F605-85ER40226, and the Vanderbilt University Research Council. 
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